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Abstract. We proceed with the construction of normalizable Dirac wave packets for 
treating chiral oscillations in the presence of an external magnetic field. Both chirality 
and helicity quantum numbers correspond to variables of fundamental importance 
in the study of chiral interactions, in particular, in the context of neutrino physics. 
In order to clarify a subtle aspect in the confront of such concepts which, for 
massive particles, represent different physical quantities, we are specifically interested 
in quantifying chiral oscillations for a fermionic Dirac- type particle (neutrino) non- 
minimally coupling with an external magnetic field B by solving the correspondent 
interacting Hamiltonian (Dirac) equation. The viability of the intermediate wave 
packet treatment becomes clear when we assume B orthogonal/parallel to the direction 
of the propagating particle. 
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Since those early days when Dirac had derived the relativistic wave equation for a 
free propagating electron [1] , several efforts have been produced in the literature to solve 
the Dirac equation with other analytical forms of interacting potentials, from central 
potential solutions [2, 3] to recent theoretical attempts to describe quark confinement 
[4, 5]. In fact, obtaining exact solutions of a generic class of Dirac wave equations 
becomes important since, for many times, the conceptual understanding of physics 
can only be brought about by such solutions. These solutions also correspond to 
valuable means for checking and improving models and numerical methods for solving 
complicated physical problems. In the context in which we intend to explore the Dirac 
formalism, we can take into account the effect of neutrino spin flipping attributed to some 
dynamic external [6] interacting processes which come from the non-minimal coupling 
of a magnetic moment with an external electromagnetic field [7] and which was formerly 
supposed to be a relevant effect in the context of the solar-neutrino puzzle by suggesting 
an explanation for the LSND anomaly [8, 9, 10] ||. 

In the correspondent theoretical framework, we can notice a more extended 
discussion involving the evolution equation of the covariant neutrino spin operator 
in the Heisenberg representation in presence of general external fields [16]. Such an 
analysis involves both electromagnetic and weak interactions with background matter 
[16, 17]. In a similar context, some developments considering the minimal coupling of 
an electrically charged particle with an external magnetic field were also performed for 
an electron described by the Dirac equation [18] and the theory of spin light of neutrino 
in matter and electromagnetic fields has been extensively studied [19]. 

We know, however, that independently of any external electromagnetic field, since 
neutrinos are detected essentially via V-A charged weak currents, the chiral oscillation 
mechanism by itself may even provide a solution to the posed anomaly. The point is 
that an erroneous correspondence is supposed to cause some confusion in the literature 
where the concept of helicity has been currently used in the place of chirality. In the 
standard model of flavor- changing interactions, neutrinos with positive chirality are 
decoupled from neutrino absorbing charged weak currents [24]. Consequently, neutrinos 
with positive chirality sterile with respect to weak interactions. By reporting to the 
formalism with Dirac wave packets [25, 26], which leads to the study of chiral oscillations 
[24] (in vacuum), we are now interested in appointing the modifications of the chirality 
dynamics which are observed when the neutrino electrodynamics is accurately discussed. 

In this context, our aim in this manuscript is to try to accommodate the Dirac 
wave packet formalism [25, 26] and a further class of static characteristics of neutrinos, 
namely, the (electro) magnetic moment which appears in a Lagrangian with non- minimal 
coupling. By following this line of reasoning, we are not only interested in solving a 
"modified" Dirac wave equation but also in constructing Dirac wave packets which can 

|| In fact, the advent of the neutrino physics [11, 12, 13] has stirred up an increasing number of 
theoretical works where the quantum oscillation mechanisms [14, 23, 21, 22, 15] have been deeply 
analyzed. Although clear experimental evidences are still missing, it is widely believed that neutrino 
mixing and oscillations are the basic tool for further understanding of neutrino phenomenology. 
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be used for describing the dynamics of chirality, a physical variable which is relevant in 
the context of the quantum oscillation phenomena. 

Despite their electric charge neutrality, neutrinos can interact with a photon 
through loop (radiative) diagrams. The Lagrangian for the interaction between a 
fermionic field i^(x) and an electromagnetic field written in terms of the field-strength 
tensor F» u (x) = d fl A u {x) - d v A»{x) is given by 

£ = ^(x) v [|iFW-<in)]#) + k (1) 

where x = (t,x), a^ v = | [7^,71/], the dual field-strength tensor T^ v {x) is given by 
T^ix) = I e M ^ A5 F xs \x)% and the coefficients /j, and d represent, respectively, the magnetic 
and the electric dipole moment which establish the neutrino electromagnetic coupling. 
One can notice that we have not discriminated the flavor/mass mixing elements in the 
above interacting Lagrangian since we are indeed interested in the physical observable 
dynamics ruled by the Hamiltonian 



^(^F^(x)-d^ u (x)) + h.c. 



H = ct ■ p + (3m — (3 

= efp + p[m-pi'£-B{x) + d'£-E{x)], (2) 

where, in terms of the Dirac matrices, a = Y? k =i a ^ = SjLi lolkk, (3 = j , and 
B(x) and E(x) are respectively the magnetic and electric fields. In fact, the Eq. (2) 
could be extended to an equivalent matrix representation with flavor and mass mixing 
elements where the diagonal (off-diagonal) elements described by fi itJ (m id ) and d itJ (m itJ ), 
where i, j are mass indices, would be called diagonal (transition) moments. In this 
context, for both Dirac and Majorana neutrinos, we could have transition amplitudes 
with non- vanishing magnetic and electric dipole moments [27, 28, 29]. Otherwise, the 
CP invariance holds the diagonal electric dipole moments null [30]. Specifically for 
Majorana neutrinos, it can be demonstrated that the diagonal magnetic and electric 
dipole moments vanish if CPT invariance is assumed [27]. 

Turning back to the simplifying example of diagonal moments and assuming CP 
and CPT invariance, we can restrict our analysis to the coupling with only an external 
magnetic field B(x) by setting d — 0. From this point, the expression for fj, can be 
found from Feynman diagrams for magnetic moment corrections [30] and turns out to 
be proportional to the neutrino mass (matrix), 

3eG 2>m e G n _ ,„_ in rn u 

u, = — = — m = — 1= — u B m u = 2.7 x 10 a B (3) 

P 8V2tt 2 4V27C 2 m N w 

where G is the Fermi constant and tun is the nucleon mass + . In particular, for 
m u PS 1 eV, the magnetic moment introduced by the above formula is exceedingly small 
to be detected or to affect astrophysical or physical processes. 

Since we are interested in constructing the dynamics ruled by the Hamiltonian 
of Eq. (2), we firstly observe that the free propagating momentum is not a conserved 

e^ uXS is the totally fourth rank antisymmetric tensor. 
+ We are using some results of the standard SU(2)l ® U(1)y electroweak theory [31]. 
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quantity, 

j f (p) = i{[H,p\) = Re(u)(f3V(X-B(x))) (4) 
In the same way, the particle velocity given by 

Ux) = i([H,x}) = (a) (5) 



comes out as a non-null value. After solving the x(t)(ct(t)) differential equation, it 
is possible to observe that, in addition to a uniform motion, the fermionic particle 
executes very rapid oscillations known as zitterbewegung [32] . By following an analogous 
procedure for the Dirac chiral operator j 5 , newly recurring to the equation of the motion, 
it is possible to have the chirality and the helicity dynamics respectively given by 

2i( lol5 [m-fiX-B(x)}) (6) 

\ ii ( 7o [(£ • V)(E • B(x)) + 2(E x B(x)) ■ p]> (7) 

where we have alternatively defined the particle helicity as the projection of the spin 
angular momentum onto the vector momentum, h = |S p (with p in place of p). From 
Eqs. (6-7) we can state that if a neutrino has an intrinsic magnetic moment and passes 
through a region filled by an magnetic field, the neutrino helicity can flip in a completely 
different way from how chiral oscillations evolve in time. In the non-interacting case, it 
is possible to verify that the time-dependent averaged value of the Dirac chiral operator 
7 5 has an oscillating behavior [24] very similar to the rapid oscillations of the position. 
The Eqs. (6-7) can be reduced to the non-interacting case by setting B(x) — so that 

i{[H,h\) = -{(axp).p) = (8) 

t([H,r']} = 2*m( 7 V>- (9) 

from which we confirm that the chiral operator 7 s is not a constant of the motion [24]. 
The effective value of Eq. (9) appears only when both positive and negative frequencies 
are taken into account to compose a Dirac wave-packet, i. e. the non-null expectation 
value of (70 7s) is revealed by the interference between Dirac equation solutions with 
opposite sign frequencies. The effective contribution due to this interference effect lead 
us to report to the Dirac wave-packet formalism in order to quantify neutrino chiral 
oscillations in the presence of an external magnetic field. 

Assuming the simplifying hypothesis of a uniform magnetic field B, the physical 
implications of the non-minimal coupling with an external magnetic field can then be 
studied by means of the eigenvalue problem expressed by the Hamiltonian equation 

H(p) ip n = E n ( P ) ip n = {a ■ p + [m - /i S • B}} <p n , (10) 



-<7 5 > = 
dt K 1 1 

and 



and 

-<7 5 > = 
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for which the explicit 4x4 matrix representation is given by 

m - fiB, -fi(B x - iB y ) p z p x - ip y 

-fj,(B x + iB y ) m + [iB z p x + ip y -p z 

p z Px-iPy -(m-fxB z ) n(B x -iB y ) 

p x + ip y -p z fi(B x + iB y ) -(m + iiB z ] 

The most general eigenvalue (E„(p)) solution of the above problem is given by 



1, 2, 



(12) 



E n (p) = ±E„{p) = ±^m 2 + p 2 + a 2 + (-l) s 2^/m 2 a 2 + (p x a) 2 , s 

where we have denoted E n=1A3j4 = ±E s=lj2 and we have set a = ji B. The complete set 
of orthonormal eigenstates <p n thus can be written in terms of the eigenfunctions U(p s ) 
with positive energy eigenvalues (+E s (p)) and the eigenfunctions V(p s ) with negative 
energy eigenvalues (—E s (p)), 

\ t 



IA- 



U{p.)= -N{ Ps ) {J-*-, J — 



V(p.) = - N( Ps ) 



s 

/A 1 

s 

1^' 



s s 

A+a-' 



1 




(13) 



where p s is the relativistic quadrimomentum, p s = (E s (p),p), N(p 3 ) is the normalization 
constant and 



Af = A 2 s (p) ± 2m\a\ - a 2 , af = 2E,(p)\a\ ± (A 2 (p) + a 2 ) 



(14) 



with 



A:(p) = £ s 2 (p)-(m 2 + p 2 ) + a 2 . (15) 

We can observe that the above spinorial solutions are free of any additional constraint, 
namely, at a given time t, they are independent functions of p and they do not represent 
chirality/helicity eigenstates. 

In order to describe the above Hamiltonian dynamics for a generic observable 0(t) 
we can firstly seek a generic plane-wave decomposition as 

exp [—i(E s (p) t — p ■ x)} U{p s ), for positive frequencies and 
exp [ i(E s (p) t — p ■ x)} V(p a ), for negative frequencies, 

so that the time-evolution of a plane-wave-packet ip(t,x) can be written as 

d?p 



(16) 



4>(t,x) 



Y,{b{Ps)U{ Ps ) exp [-iE,(p)t] 



+ d*(p s )V(p s ) exp[+iE.(p)t]}exp[ip-x], (17) 

with p s = (E s ,—p). Meanwhile, the Eq. (17) requires some extensive mathematical 
manipulations for explicitly constructing the dynamics of an operator 0(t) like 



0(t) 



d 3 xip\t,x)(D ip(t,x). 



If, however, the quoted observables like the chirality j 5 , the helicity h or even the spin 
projection onto B commuted with the Hamiltonian H, we could reconfigure the above 
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solutions to simpler ones. To illustrate this point, let us limit our analysis to very- 
restrictive spatial configurations of B so that, as a first attempt, we can calculate the 
observable expectation values which appear in Eq. (6). Let us then assume that the 
magnetic field B is either orthogonal or parallel to the momentum p. For both of these 
cases the spinor eigenstates can then be decomposed into orthonormal bi-spinors as 



Ufa) = N + ( Ps ) 



X + (Ps) 



(19) 



and 



V(p.) = N-( Ps ) 



(20) 



V (Ps) 
X~(Ps) 

Eventually, in order to simplify some subsequent calculations involving chiral 
oscillations, we could set ipf 2 and xf 2 as eigenstates of the spin projection operator 
cr ■ B, i. e. beside of being energy eigenstates, the general solutions U(p s ) and V(p s ) 
would become eigenstates of the operator £ ■ B and, equivalently, of £ • a. 
Now the Eq.(10) can be decomposed into a pair of coupled equations like 

(±E,-m + <r-a)(pf = ±<r-pxt, 



m 



a )xf = ±<T-p<p* 



(21) 



where we have suppressed the p s dependence. By introducing the commuting relation 
[cr • p, cr ■ B] = which is derived when p x B = 0, the eigenspinor representation can 
be reduced to 



and 



W(p.) = 



V(p.) 









/ 2E S 


cr p 

L E s +m s 





l E a + m s 



IE, 



cr p 



Xs 



with m s — m — (-l) s |a| and the energy eigenvalues 

±E S = ±^ + mJ. 
In this case, the closure relations can be constructed in terms of 



(22) 



(23) 



(24) 



u{ Ps )®u\ Ps ) lo 



8 = 1,2 



E v(p.) ® v\ Ps ) lo 




2E S 
IE, 



1 - (-1) S E • a 



1 - (-1) S S • a 



(25) 



Analogously, by introducing the anti-commuting relation {cr • p, cr ■ B} when p ■ B = 0, 
the eigenspinor representation can be reduced to 



K(Ps) 



le + m 






1 2e 


<T p 

- eo+m 





(26) 
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V(p.) 



je + m 


r 




1 2 £o 




x; . 



with e 



^/p 2 + m 2 and the energy eigenvalues 
±£. = ± [e + (-l)'|a|] . 
In this case, the closure relations can be constructed in terms of 

"1 - (-l)' 7o S • o' 



(27) 



(28) 



s=l,2 s= i,2 



- E V(p.) <8> V t (p 3 )7o 



~7mPu + m 
2e 



E 



1 - (-l) s 7oS • a 



(29) 



where p = (e ,p). 

Since we can set t/?+ 2 = x7, 2 as the components of an orthonormal basis, we can 
immediately deduce the orthogonality relations 

UH P s)K(Pr) =V\p 3 )V( Pr ) = 5 S „ 

U\ Ps ) 7o V(p r ) = V f ( Ps ) 7o W(p r ) = (30) 

which are valid for both of the above cases. 

Finally, the calculation of the expectation value of j 5 (t) is substantially simplified 
when we substitute the above closure relations into the wave-packet expression of 
Eq. (17). To clarify this point, we suppose the initial condition over ip(t,x) can be 
set in terms of the Fourier transform of the weight function 

<p(p-Pi)w = E {b(p 8 )K(p s ) + d*(ps)V(p s )} (31) 



so that 



/d l p 
j2nyM p ~ p ^ ex P [ i P- x ) w 



(32) 



where w is some fixed normalized spinor. By using the orthonormality properties of 
Eq. (30), we find [25] 

b(p 3 ) = ip{p- Pi )U\p s )w, 

d*(p s ) = ip(p- Pi )VHPs)w. (33) 

For any initial state ip(Q,x) given by Eq. (32), the negative frequency solution coefficient 
d*(p a ) necessarily provides a non-null contribution to the time-evolving wave-packet. 
This obliges us to take the complete set of Dirac equation solutions to construct a 
complete and correct wave-packet solution. Only if we consider the initial spinor w 
being a positive energy (E 3 (p)) and momentum p eigenstate, the contribution due to the 
negative frequency solutions d*(p~ s ) will become null and we will have a simple expression 
for the time-evolution of any physical observable. By substituting the closure relations 
of Eqs. (25) and (29) into the time-evolution expression for the above wave-packet, the 
Eq.(17) can thus be rewritten as 



(2tt> 



-¥>(p-Pi)exp [ip-x] 



3=1,2 



l-i-iy'E-d 

2 



w 



(34) 
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for the first case where E s is given by Eq. (24) and H s = a ■ p + %m s , or as 



i/j(t,x)-- 



J T^y v[p Pi) exp [ip ' x] 5 { H [Es t] ~ sin [E ° t] \ ( 1 ' ( ' 1) 2 7oE& ) } w ( 35 ) 

for the second case where e is given by Eq. (28) and H — a • p + j m. 

Once we have assumed the neutrino electroweak interactions at the source and 
detector are (left) chiral (-?/>7 M (i — y')ipWA only the component with negative chirality 



contributes to the propagation. It was already demonstrated that, in vacuum, chiral 
oscillations can introduce very small modifications to the neutrino conversion formula 
[24, 20]. The probability of a neutrino produced as a negative chiral eigenstate be 
detected after a time t can be summarized by the expression 

P(v a ,L - u^-t) = J d 3 xip\t,x)^^ij(t, x ) = l - (1 - ( 7B )(t)) . (36) 

From this integral, it is readily seen that an initial -1 chiral mass-eigenstate will evolve 
with time changing its chirality. By assuming the fermionic particle is created at time 
t = as a -1 chiral eigenstate ij 5 w = -w), in case of [cr • p, a • B\ — (B parallel to 
p), we could write 



i-(-i) s E-q 

2 



™ f £ { h c «s 2 [E. t] + sin [2 E s t] + sin 2 [E. t] 

3 = 1,2 ^ 

= (-D / - pj E { [ £ . *i + ^ ™ a [ £ . <]] ( M i' s ° 

where we have used the wave packet expression of Eq. (34) and, in the second passage, 
we have observed that 

1^7,. w = -1, w^[H s , j 5 ]w = and # s 7 5 H s = p 2 -m 2 s . (38) 

The above expression can be reduced to a simpler one in the non-interacting case [24]. 
Due to a residual interaction with the external magnetic field B we could also observe 
chiral oscillations in the ultra-relativistic limit. However, from the phenomeno logical 

2 

point of view, the coefficient of the oscillating term goes with which makes chiral 
oscillations become not relevant for ultra-relativistic neutrinos [20, 15]. As a toy model 
illustration, by assuming a highly peaked momentum distribution centered around a non- 
relativistic momentum p t <C m s , where the wave packet effects are practically ignored, 
the chiral conversion formula can be written as 

P{v a , L — > Va,ii] t) ^ (l — cos [2mt] cos [2 \ a\t] — sin [2mt] sin [2 \ a\ t]w^H ■ awj (39) 

where all the oscillating terms come from the interference between positive and negative 
frequency solutions which compose the wave packets. Turning back to the case 
where {er ■ p, a ■ B} — (B orthogonal to p), we could have a phenomenologically 
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more interesting result. By following a similar procedure with the mathematical 
manipulations, we could write 

(75>(*) = / J^yytP Pi) w *{l* cos [^i *] cos i E * *] + sin i E i *] sin *] + 

f [[#„, 7b ] sin + E 2 ] - {H , 7b } 7o S • a sin - £ 2 ]]} w 
= (-1) / (07^ 2 (P - Pi){cos t] cos [E 2 1] + S i n [ El t] sin [£ 2 *]} 

= (-1) /?0i^-^){f cos[2|a|t] + ^cos[2e t]}, (40) 

where we have used the wave packet expression of Eq. (35) and, in addition to w'j s w = 
-1, we have also observed that {H , 7 5 } = 27 5 S-p and, subsequently, w^Yi-fry Yi-aw = 0. 
Now, in addition to the non-interacting oscillating term ^- cos [2 e i] , which comes from 
the interference between positive and negative frequency solutions of the Dirac equation, 
we have an extra term which comes from the interference between equal sign frequencies 
and, for very large time scales, can substantially change the oscillating results. In this 
case, it is interesting to observe that the ultra-relativistic limit of Eq.(40) leads to the 
following expressions for the chiral conversion formulas, 

P{v a ,L -> t) w - (1 - cos [2 \a\ t}) (41) 

and 

P{v a ,L -> t) w X - (1 + cos [2 |a| *]) (42) 

which, differently from chiral oscillations in vacuum, can be phenomeno logically relevant. 
Obviously, we are reproducing the consolidated results already attributed to neutrino 
spin-flipping [30] where, by taking the ultra-relativistic limit, the chirality quantum 
number can be approximated by the helicity quantum number, but now it was accurately 
derived from the complete formalism with Dirac spinors. We still remark that, in the 
standard treatment of vacuum neutrino oscillations, the use of scalar mass-eigenstate 
wave packets made up exclusively of positive frequency plane-wave solutions is implicitly 
assumed. A satisfactory description of fermionic (spin one-half) particles requires the use 
of the Dirac equation as evolution equation for the mass-eigenstates despite the standard 
oscillation formula giving the correct result when it is properly interpreted. We have 
observed in Eqs. (40) and (37) that the spinorial form and the interference between 
positive and negative frequency components of the mass-eigenstate wave packets can 
introduce small modifications to the standard conversion formulas when they concern 
with non-relativistic neutrinos [15]. In the next step, by following the Dirac wave packet 
prescription, we intend to study the coupling of chiral oscillations with flavor conversion 
effects for neutrinos non-minimally interacting with an external magnetic field and, 
subsequently, verify the possibility of some phenomenological implications. 

To summarize, we would have been dishonest if we had ignored the complete 
analysis of the general case comprised by Eqs. (10-13) where we had not yet assumed 
an arbitrary (simplified) spatial configuration for the magnetic field. Meanwhile, such 
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a general case leads to the formal connection between quantum oscillation phenomena 
and a very different field. It concerns with the curious fact that those complete (general) 
expressions for propagating wave packets do not satisfy the standard dispersion relations 
like E 2 = m 2 + p 2 excepting by the two particular cases where E a (p) 2 = m 2 + p 2 for 
p x B = or el = m 2 +p 2 for p -B = 0. By principle, it could represent an inconvenient 
obstacle forbidding the extension of these restrictive cases to a general one. However, 
we believe that it can also represent a starting point in discussing dispersion relations 
which may be incorporated into frameworks encoding the breakdown (or the validity) 
of Lorentz invariance. 
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